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This paper presents a novel self-tuning fractional order PID (FOPID) control 
based on optimal Model Reference Adaptive Control (MRAC). The proposed 
control technique has subjected to a third order system case study (power 
system load frequency control). The model reference describes the 
requirements of designer. It can be first or second order system. The 
parameters of MRAC have obtained using the harmony search (HS) 
optimization technique to achieve the optimal performance. Sometimes, the 
tuning of the five parameters of FOPID control online at same moment 
consumes more calculation time and more processing. So, this study 
proposes three methods for self-tuning FOPID control. The first method has 
been implemented to tune the two integral and derivative parameters only 
and the rest of parameters are fixed. The second method has been designed to 
adjust the proportional, integral derivative parameters while the other 
fractional parameters are constant. The last method has developed to adjust 
the five parameters of FOPID control simultaneously. The simulation results 
illustrate that the third method of self-tuning FOPID control can 
accommodate the sudden disturbance compared to other techniques. Also, it 
can absorb the system uncertainty better than the other control techniques. 
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1. INTRODUCTION 

The PID control has been applied to process control in most of engineering applications for decades 
[1], The PID control has simple structure and linear behavior. Also, it gives acceptable performance for 
several industrial applications [2], There are several methods to select the proper values for PID controller 
parameters [3], The traditional methods for selecting these parameters such as try and error and Ziegler- 
Nichols which were became inappropriate to achieve a good performance [4], So, the researchers have 
tended to use alternative methods such as optimization techniques (Genetic Algorithm (GA), Particle Swann 
Optimization (PSO), Ant Colony Optimization (AGO) and Hannony Search (HS)) which are trying to reach 
the optimal solution for controller parameters [5], Still, the behavior of PID control is linear and cannot deal 
with the high disturbance and high nonlinearity of complicated systems [6], [7], 

The fractional order PID (FOPID) control has been widely used in control engineering in recent 
decades [8], The FOPID considers the nonlinear copy of PID control where two more parameters (the 
fractional integral and derivative) added to the PID control parameters [9], Hence, the task of designer 
selecting the proper values for the five parameters of the FOPID control [10]. The FOPID control can solve 
the nonlinearity problem but it cannot deal with the sudden disturbance due to its parameters which 
still fixed [11], 
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Different techniques have been proposed to tune the five parameters of FOPID control online but 
most of these techniques are based on the fuzzy logic control [12], [13]. The fuzzy logic control can solve the 
uncertainty problem and sudden disturbance but its design depends on the experience which sometimes is not 
available for some systems [ 14]-[ 16]. This study presents a novel technique to tune the FOPID control 
parameters online based on optimal model reference adaptive control (MRAC). It is known that the MRAC is 
high ranking adaptive control where it forces the overall system to follow the behavior of preselected model 
reference [17]. The preselected model can be first or second order system according to the point of view the 
designer and complicated degree of the system [18]. 

The task of model reference adaptive control is adjusting the FOPID control parameters online. The 
model reference contains the desired performance which can satisfy the designer. Moreover, to guarantee 
high performance the parameters of model reference optimized using the harmony search (HS) optimization 
technique according to a certain cost function. 

The proposed technique will be subjected to a third order system as case study (power system load 
frequency control). Also, the proposed technique has been implemented with different methods. The first 
method has been designed to the MRAC will tune the two integral and derivative parameters only and the 
rest of parameters are fixed. The second method has been developed to adjust the proportional, integral 
derivative parameters while the other fractional parameters are constant. The last method fabricated to adjust 
the five parameters of FOPID control at the same time. The paper has organized as follows, firstly, the 
system model is presented. Secondly, the proposed control techniques are demonstrated. Thirdly, the 
simulation results are illustrated. Finally the conclusion is discussed. 


2. THIRD ORDER CASE STUDY MODEL 

Several nonlinear models for large power systems has to be established but all are complicated 
models. So, the linearized model has been used in our work. The power system represented by an equivalent 
turbine, governor and generator system [19], [20]. The transfer function of a single area power system model 
in a closed loop form as shown in Figure 1. 



Figure 1. Block diagram of a single area power system 


Where is u the control action, d is the disturbance vector (APd), Tp is the plant model time constant, 
Tt is the turbine time constant, Tg is the governor time constant, Kp is the plant gain, R is the speed 
regulation due to governor action, xl is the change in system frequency, x2 is the incremental changes in 
generator output and x3 is the governor valve position. Table 1 present the values of system parameters that 
used in this work. 


Table 1. System Parameters 


Parameters 

Value 

Kp 

120 pu 

Tp 

20 s 

Tt 

0.3 s 

Tg 

0.08 s 

R 

2.4 Hz/pu.MW 


3. SELF-TUNING FOPID CONTROL 

The Model Reference Adaptive Control (MRAC) is high-ranking adaptive controller [2], [5]. It may 
be regarded as an adaptive servo system in which the desired perfonnance is expressed in tenns of a 


A novel self-tuning fractional order PID control based on optimal model... (Moliamed. A. Shamseldin) 











232 □ 


ISSN: 2088-8694 


reference model. In this work the FOPID control parameters will be adjusted on-line using the model 
reference technique as shown in Figure 2. 



Figure 2. The overall system with self-tuning FOPID based on model reference technique 


To adjust parameters using MIT rule which use the following loss function 




ID 


j, that is, 


To make j small, it is reasonable to change the parameters in the direction of the negative gradient of 


d6 dj dei 

*<■ ~ ^ ae ~ ^ e,n 


ae 


(2) 


where y stand for the adaptation gain while 0 is the centroid vector of the output membership function. The 
transfer function of FOPID control can be described as follows. 


^ = k p + k i ± + k d st l 

e(s) v 1 s A u 


(3) 


e = u c — y (4) 

Assume that the plant can be simplified to a first order system as obvious in the following (5). 


yQ) _ k 

u(s) Ts+1 


(5) 


where k and T are unknown parameters. Also, assume that the model reference takes a form first order 
system as the following relationship. 
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Where k m and T m are selected by designer.From (4-6) can conclude that 
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3.1. Adaptation law of k p parameter 

This sub-section shows the steps of design of the adaptation law for proportional gain parameter 
(k p ). By deriving (10) respect to the proportional gain ( k p ) to obtain the following relationship. 
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From (12) and (14); 
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To achieve the desired performance, the following condition must be hold. 
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dk p _ _ e m .e 

dt ~ ^ 1 'T m s+l 

(19) 
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II 

£ 

(20) 

k v )new = J^T dt + k p ( 0) 

(21) 


Where k p ( 0) is the initial value of proportional gain k p . 


3.2. Adaptation law of k t parameter 

This sub-section shows the steps of design of the adaptation law for integral gain parameter (fc;). By 
deriving (8) respect to the integral gain (fc;) to obtain the following relationship. 
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To achieve the desired performance, the condition must be hold in (16). 

de m _ 1 k 2 e 
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From the MIT rule can obtain the following relationship 
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Where fcj( 0) is the initial value of proportional gain k t . 
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3.3. Adaptation Law of k d Parameter 

This sub-section illustrates the steps of design of the adaptation law for derivative gain parameter 
( k d ). By deriving (10) respect to the derivative gain ( k d ) to obtain the following relationship. 
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Also, from (33) and (35) 
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Where k d ( 0) is the initial value of derivative gain k d . 
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3.4. Adaptation Law of A Parameter 

This sub-section illustrates the steps of design of the adaptation law for fractional integral gain 
parameter (2). By deriving (10) respect to the fractional integral gain (A) to obtain the following relationship. 
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de m _ kkiln(s) 
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Also, from (41) and (43); 
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Where 1(0) is the initial value of fractional integral gain X. 
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3.5. Adaptation law of/r parameter 

This sub-section demonstrates the steps of design of the adaptation law for fractional derivative gain 
parameter (n). By deriving the (8) respect to the fractional derivative gain (/ 1 ) to obtain the 
following relationship. 
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Also, from (51) and (53); 
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_ _ £mf 

dt ^ 5 'r m s+i 

Ts = y.k 2 .k d (0).s^ o \ln(is) = y 3 ./c d (0)./n(s) 

M) new = J^r dt + M(0) 

Where fi(0) is the initial value of fractional integral gain p. 

3.6. Harmony search (HS) optimization 

The harmony search has been proposed in 2001. It is used lately in several engineering applications 
to obtain the optimal values of control parameters. In this study, the HS will be used to obtain the optimal 
parameters values of model reference adaptive control. The most effective parameter of MRAC is the 
adaptation gain for each adaptation law. Here, the harmony search tuning system will adjust the parameters 
of model reference adaptive system (7i, 72 - 73-74 an d 7s) according to the objective function as 
shown in (60). 

f (1 -e-P)(M v +e ss )+e~P(t s -t r ) ^ 0 ) 


(57) 

(58) 

(59) 


Where M p the overshoot of system response is, e ss is the steady state error, t s is the settling time and t r is the 
rise time. Also, this objective function is able to compromise the designer demand by the weighting 
parameter value ((3). The parameter is set larger than 0.7 to reduce over shoot and steady-state error. If this 
parameter is set smaller than 0.7 the rise time and settling time will be reduced. 

The initial population of Hannony Memory (HM) is produced randomly. HM contains Harmony 
Memory Solution (HMS) vectors. The HM is filled with HMS vectors as follows: 


Yi 

(i.i) 

72 (1,2) 

Ys (1,3) 

Y* (1,4) 

7s 

Yi 

(2,1) 

Y2( 2,2) 

^3(2,3) 

Y 4 (2,4) 

7s 


HM = 


(61) 


lYl (HMS, 1)^2 (HMS,2)Y’S(HMS,S)Ya (HMS,4)7S( HM s, 5 )J 

Table 2 demonstrates the obtained values of MRAC parameters after the offline tuning using the 
harmony search tuning system. 


_ Table 2. MRAC Parameters _ 

mrac Yi _ y-/ y-a ft _ 7a 

Parameters 0.234 0.678 0.456 0.568 0.7043 


Sometimes, some systems don’t need adjust the five parameters of FOPID control online at same 
moment which will save the calculation time and make the overall system more ready to deal with the sudden 
disturbance. So, this study proposes three methods for self-tuning FOPID control. The first method 
considers k p , k t and k d are constants while X and p are varying. The second method considers k p , k t and 
k d are varying while X and p are constants. The third method tunes the five parameters of FOPID control 
online simultaneously. 


4. SIMULATION RESULTS 

This section presents the simulation results of the proposed different types of self-tuning FOPID 
control algorithms based on model reference. The first method, the adaptive mechanism will adjust the 
fractional integral and derivative parameters only while the other parameters are fixed. The second method, 
the fractional integral and derivative parameters are fixed while the other three parameters will be tuned 
using the adaptive mechanism. The third method, the five parameters of FOPID control will be adjusted 
simultaneously. Figure 3 demonstrates the Simulink diagram of the overall system with self-tuning FOPID 
based on optimal model reference technique. 
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Figure 3. The simulink diagram of the overall system with self-tuning FOPID based on optimal model 

reference technique 


The results are divided into two cases. The normal case where applied two types of the disturbers. 
These disturbers namely APdl as shown in Figure 4 and APd2 as illustrated in Figure 5. The parametric 
uncertainty case where to change the operating point to test the powerful of the proposed method against 
their counterparts 
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Figure 4. Power demand variation (APdl) 


Figure 5. Power demand variation (APd2) 


4.1. Normal case 

Two experiments were done at normal operating condition. The dynamic response of the system 
frequency AF with the system driven by each type of the proposed controllers. In experiment number one 
where APdl is applied to the system. Figure 6 shows the time response of the system. Also, the output signal 
of each controller is shown in Figure 7. In experiment number two, Figure 8 demonstrates the time response 
of the system when APd2 is applied to the system. Figure 9 presents the controller’s outputs signals. 
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Figure 8. System dynamic responses at APd2 


0.3- 



Figure 9. The controller's output signals at APd2 


It is clearly shown that The proposed third method of self-tuning FOPID controller gives battery 
performance with a smaller settling time and acceptable undershoot but with a relatively higher effort as 
shown in the response of its control output response u. also, the proposed third type of self-tuning FOPID 
controller recover to zero steady state error after a smaller time delay from the application of the two 
disturbance types. 

4.2. Parametric uncertainty case 

In this subsection, the parametric uncertainties of the power system have to be considered. 
According to load variation and power system configuration, the operation points of the system will be 
changed randomly during a daily cycle. The system parametric uncertainties are obtained by changing 
parameters by 50% from their normal values according to Table 3 and under the power demand variation 
APd2. In this experiment, the power system responses at APd2, including the effect of the parametric 
uncertainties is presented in Figure 10 and the controllers outputs are shown in Figure 11. 

According to these results, it is clear that the proposed third method of self-tuning FOPID control 
shows the best response compared to other controllers. Also, the results approve the effectiveness and the 
ability of the proposed controller against the parametric uncertainties. Moreover, the undershoot at 1 second 
and overshoot at 5 second of third method of self-tuning FOPID control are significant small compared to 
other methods of self-tuning FOPID control. 
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Table 3. Parametric Uncertainties 


Parameters 

Down 
(0-1) sec 

Normal 
(1-5) sec 

Up 

(5-7)sec 

k p 

60 

120 

180 

T P 

10 

20 

30 

R 

1.2 

2.4 

3.6 

7> 

0.15 

0.3 

0.45 

To 

0.04 

0.08 

0.12 




Figure 10. System dynamic responses at APd2 and Figure 11. The controller's output signals at APd2 and 
parametric uncertainties parametric uncertainties 


5. CONCLUSION 

A novel technique has been developed to tune the fractional order PID (FOPID) control online 
based on optimal Model Reference Adaptive Control (MRAC) for power system load frequency control as a 
case study. This work investigates the robustness of the proposed technique by applied different types of 
disturbances and parametric uncertainty. There are three methods of self-tuning FOPID control have been 
implemented. The first method, the three integer parameters (proportional, integral and derivative gains) are 
constant while the fractional integral and derivative gains are variable. The second method, the fractional 
integral and derivative gains are constant while the three integer parameters (proportional, integral and 
derivative gains) are variable. The third method, the all five parameters are variable and change continuously 
online. The simulation results demonstrate that the third method of self-tuning FOPID control can absorb the 
sudden disturbance faster than other techniques. Also, it can suppress the system uncertainty in a small time 
compared to other control techniques. 
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